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EXECUTIVE  SUMMARY 


INTRODUCTION 

The  Adaptive  Optics  task  was  funded  under  the  IR  Program.  The  purpose 
of  the  study  was  twofold: 

a.  Determine  the  importance  of  the  air/sea  interface  in  the  degradation 
of  laser  beam  irradiance 

b.  Examine  the  possibility  of  using  adaptive  optical  techniques  to  com- 
pensate for  the  interface 

BACKGROUND 

When  a laser  beam  is  directed  from  an  airborne  platform  to  a target  under 
the  water,  the  irradiance  is  diminished  by  the  initial  beam  divergence,  atmos- 
pheric scattering  and  absorption,  the  air/sea  interface,  and  bulk  scattering  in 
the  ocean.  By  using  the  small  angle  scattering  approximation  of  electron 
scattering  theory  and  a modified  Cox  and  Munk  probability  distribution  for  sea 
slopes,  an  expression  is  derived  which  predicts  the  laser  beam  radius  as  a 
function  of  distance  (z)  into  the  water.  A comparison  of  the  terms  in  the 
resultant  expression  for  beam  radius  shows  that  the  bulk  scattering  term 
dominates  over  the  surface  spreading  term  for  z>10  meters  in  very  good  water 
(s~.  05m  S and  high  wind  speeds  (~10m/s). 

It  was  decided  to  make  a caretui  determination  of  the  capillary  wave  slope 
distribution  in  the  laboratory  because  the  Cox  and  Munk  values  are  smaller  than 
the  theoretical  values.  Additionally,  if  the  slope  distribution  was  determined 
to  be  much  larger,  then  the  conclusion  about  the  relative  importance  of  the 
interface  would  be  changed  accordingly. 

CONCLUSIONS 

The  following  conclusions  can  be  drawn  from  this  study: 

a.  The  beam  spreading  due  to  the  roughness  of  the  air-sea  interface  is 
small  compared  with  the  spreading  due  to  bulk  scattering.  Therefore, 
compensating  for  the  interface  would  be  of  little  benefit. 
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b.  Using  a small-diameter  laser  beam  substantially  increases  the  irradi- 
ance  on  target  for  depths  of  50  meters  and  less. 

c.  The  mean  heights  for  capillary  waves  are  in  the  range  0. 06  to  8. 8 milli- 
meters for  wind  speeds  2. 2 to  7. 6 meters  per  second.  These  heights 
are  consistent  with  theoretically  allowed  values  for  capillary  waves. 

d.  The  spectral  content  of  the  wind-driven  capillary  waves  is  in  the  range 

0 to  400  Hz  with  most  of  the  intensity  below  100  Hz.  The  high-frequency 
portions  of  the  spectral  curves  grow  at  the  expense  of  the  low-frequency 
portions  as  the  vind  speed  increases. 

e.  The  wavetank  data  on  wave  slopes  yield  higher  values  for  full-width 
beam  divergence  (~5  degrees)  than  open-ocean  data  (—2  degrees  and 
—4  degrees)  at  wind  speeds  in  the  range  of  5 to  8 meters  per  second. 

f.  Fewer  wave  slopes  of  high  angle  are  measured  in  these  experiments 
than  are  predicted  by  theory. 

g.  The  onset  wind  velocity  for  wave  generation  is  found  to  be  ~3  meters 
per  second. 

h.  The  shape  of  the  wind-driven  waves  is  affected  by  the  wind. 

i.  The  wave  height  measurements  furnish  experimental  support  to  the 
theoretical  model  of  Crapper. 

RECOMMENDATIONS 

1.  It  is  not  worthwhile  pursuing  adaptive  optical  methods  to  compensate 
for  ocean  roughness. 

2.  For  the  purposes  of  air-to-underwater  communication,  bathymetry, 
and  designation,  higher  rms  wave  slopes  can  be  used  than  those  given 
by  Cox  and  Munk. 
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SECTION  1 
INTRODUCTION 


This  study  was  initiated  to  establish  the  feasibility  of  using  adaptive  optic 
techniques  to  compensate  for  the  distortion  of  liglu  (laser)  beams  caused  by 
traversing  the  air-sea  interface.  Before  designing  a specific  system  which 
could  sense  the  ocean  surface  waveforms  and  tailor  the  laser  beam  wave  front 
to  compensate  for  the  interface,  it  was  decided  to  model  the  interface,  the 
volume  scattering,  and  the  incident  beam  characteristics  to  establish  the  rel- 
ative effect  of  the  interface. 


Several  wavetank  experiments  on  wind -driven  capillary  waves  were  planned 
and  performed.  The  purpose  of  these  experiments  was  to  measure  the  capillary 
wave  heights,  slopes,  and  frequencies.  These  physical  parameters  are  impor- 
tant to  the  adaptive  optics  approach. 


The  hydrodynamics  of  waves  in  water  has  been  worked  out  by  various 
authors.  A particularly  readable  approach  has  been  published  by  Clark  (Refer- 
ence 1).  A result  of  this  calculation  is  that  the  speed  (U)  of  propagation  for  a 
gravity  wave  is : 


where: 


tanh  kz 


CO 

k 

g 

z 


the 

the 

the 

the 


frequency, 
wave  vector  (^f-) , 
gravitational  constant,  and 
depth . 


(1) 
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Where  A<  < z,  equation  (1)  becomes 
1/2 

u*(k)  • <2> 

i /s 

In  accordance  with  equation  (2),  the  velocity  of  gravity  waves  increases  as  X 1 
power.  Stokes  (Reference  2)  showed  that  the  maximum  height,  trough  to  crest, 
for  gravity  waves  is  0. 14  wavelengths.  If  the  surface  tension  (S)  is  included  in 
the  calculation  for  wave  velocity,  we  find  that: 


r*  “i  l/2  l/2 

u = [(f) + (v)]  (taDh  k2) 

where  p is  the  density.  Agair,  where  A<<  z,  equation  (3)  becomes: 


(3) 


(4) 


A plot  of  equation  (4)  is  shown  in  Figure  1.  From  the  curves  in  Figure  1, 
it  can  be  seen  that  the  minimum  velocity  fcr  capillary  waves  is  approximately 
23  cm/sec;  the  minimum  frequency  is  about  13  Hz;  and  the  maximum  wavelength 
is  1.7  cm.  Also,  it  is  seen  that  the  wave  velocity  increases  as  the  capillary 
wavelength  decreases  (anomalous  dispersion)  in  contrast  to  the  gravity  wave 
velocity.  One  additional  feature  of  the  frequency  curve  shown  in  Figure  1 is 
the  rapid  increase  in  frequency  with  a decrease  in  the  capillary  wavelength.  In 
a theoretical  paper  by  Crapper  (Reference  3),  it  is  shown  that  the  wave  of 
greatest  height  occurs  when  the  trough  to  crest  distance  is  0.73  X.  This  num- 
ber is  approximately  five  times  greater  than  the  comparable  one  for  gravity 
waves.  It  is  also  shown  in  this  paper  that  the  velocity  of  capillary  waves  is 

d. 

moderately  diminished  at  higher  amplitudes  (the  factor  being  0.8  for  j = 0.73). 

The  slope  of  a wave  at  any  point  is  important  to  this  study  because  of  the 
relationship  between  slope  angle  and  beam  deflection.  The  change  of  slope  over 
the  laser  beam  area  causes  the  beam  to  spread.  Therefore,  it  is  essential  to 
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understand  wave  slope  statistics.  Of  the  many  studies  on  wave  slope  statistics 
in  the  open  ocean,  the  Cox  and  Munk  study  (Reference  4)  is  the  most  compre- 
hensive. One  of  the  results  of  the  Cox  and  Munk  study  is  that  the  wave  slopes 
may  be  described  by  a Gram  Charlier  distribution.  This  distribution  is  the 
product  of  a two-dimensional  Gaussian  distribution  and  an  infinite  series  of 
Hermite  polynomials.  To  a very  good  approximation,  the  Gram  Charlier  dis- 
tribution may  be  replaced  by  a Gaussian  distribution  as  follows: 


1 

exp  - 

2 

<T 

/ 1 + . 

1 N 

2 ir  *c  ffu 

T 

u 

_ 

2 2 

The  mean  square  slope  components  are  <?u  and  ac  , the  upwind  and  crosswind, 
respectively.  Cox  and  Munk  measured  these  slopes  to  be: 

«■  2 = 0.003  + 1.92  x 10“3  W + 0.002  ffit 

c — 

and 

"u2  = 0.000  + 3.16  x 10"3  W + 0.004  (7) 


where  W is  the  wind  speed  in  meters  per  second  (m/sec).  For  the  one- 
dimensional case,  equation  (5)  can  be  reduced  to; 


o 


where  °Q  = <7U . The  probability  distribution  in  equation  (8)  is  normalized  to 
one  (1);  that  is, 


d <r 


= 1 . 


(9) 


Two  features  of  the  above  distribution  are  that  it  is  very  narrow  at  low  wind 
speeds,  and  the  average  wave  slope  is  zero.  Neither  feature  is  exactly  true, 
but  they  are  reasonable  approximations.  The  wave  slope  distribution  can  now 
be  used  in  the  calculation  for  beam  steering  and  beam  spreading. 
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SECTION  2 
SMALL  BEAM  CASE 


Consider  a laser  beam,  small  in  diameter  compared  to  the  length  of  the 
refracting  ocean  wave,  incident  on  a wave  of  slope  a.  As  depicted  in  Figure  2, 
the  transmitted  beam  has  deviated  from  its  original  direction  by  an  amount  A. 
The  relationship  between  a and  A is  determined  by  Snell's  Law  as  follows: 

n sin  it  = n'  sin  0'  (10) 

and,  in  the  small  angle  approximation, 

n a » n'  0'.  (11) 

From  Figure  2, 


or 


For  the  case  of  seawater, 

A - ^ . (14) 

From  equation  (14)  it  can  be  seen  that  the  beam  deflection  (A)  is  approximately 
one-fourth  of  the  wave  slope  (<r).  This  result  will  be  used  in  the  one-dimensional 
wave  slope  distribution,  equation  (8),  to  determine  beam  deflection  probability. 

Substituting  equation  (14)  into  equation  (8)  and  normalizing,  we  find  that: 
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Figure  2,  Beam  Deflection  Geometry 
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Equation  (15)  shows  the  probability  that  the  beam  is  deflected  by  an  amount  A. 

In  order  to  determine  the  angular  limits  that  contain  the  beam  for  a given  per- 
centage of  the  time,  equation  (15)  is  integrated.  The  results  of  this  integration 
are  given  in  Table  1.  The  significance  of  the  results  in  Table  1 is  that  the  beam 
deflection  is  related  to  the  root  mean  square  (rms)  wave  slope  aQ  by  a probabil- 
ity. For  example,  the  probability  is  0.  8 that  A <0.23  a , and  the  probability 
is  0.5  that  A <0.12  o . For  this  calculation,  A is  the  half-angle  value.  The 
results  of  calculating  the  probability  that  the  beam  deflection  is  within  defined 
angular  limits  at  a given  wind  speed  are  listed  in  Table  2.  For  example,  for  a 
4-meter-per-second  wind  speed,  the  probability  is  0.50  that  the  beam  is  deflected 
within  13  milliradians  and  0.80  that  the  beam  is  deflected  within  26  milliradians. 
The  results  in  Table  2 for  the  small  beam  case  will  be  compared  with  the  beam- 
spreading results  for  the  large  beam  case. 


TABLE  1.  INTEGRATED  PROBABILITY  VERSUS  A/a 

o 


INTEGRATED 

PROBABILITY 

0 

0.2 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

0.96 

A/<7 

O 

0 

0.046 

0.094 

0.12 

0.15 

0.184 

0.23 

0.29 

0.36 

TABLE  2.  BEAM  DEFLECTION  VERSUS  WIND  SPEED 


WIND  SPEED 
(m/sec) 

°o 

(radians) 

A (P  = 0.5) 
(radians) 

A(P  = 0.8) 
(radians) 

4 

0.11 

0.013 

0.026 

8 

0.16 

0.019 

0.037 

12 

0.19 

0.023 

0.045 
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SECTION  3 
LARGE  BEAM  CASE 


For  the  case  in  which  the  diameter  of  the  laser  beam  is  large  compared 
with  the  ocean  wavelength  encountered,  the  beam  is  spread  by  the  air/sea  inter- 
face as  opposed  to  being  deflected.  The  large  beam  case  is  the  integral  of  the 
small  beam  case  over  the  beam  diameter.  We  assume  that  the  input  beam  f 

o 

(0,  r)  is  distributed  in  the  following  manner: 


fo  <*.  r>  - ttt 

o o 


exp 


-Bt)’ 


2-i 


R 


(16) 


where 


and 


R0  = the  radius  at  which  the  beam  decays  to  - times  its  value  at  r = o 

G 

= the  angle  at  which  the  beam  decays  to  ^ times  its  value  at  0 - o. 


The  output  beam  distribution  function  f^  (O',  r),  which  describes  how  the  beam 
is  distributed  under  the  sea  surface,  is  obtained  by  integrating  the  product  of 
the  input  beam  distribution  function  and  a modified  probability  distribution  for 
sea  slopes  over  all  input  angles,  as  follows: 


(O',  r)  = j fQ  (0,  r)  P(0\  0)  d0  . 


(17) 


In  equation  (17),  6 and  O'  are  measured  from  the  vertical,  and  P ( 0 ',  0 ) is  the 
probability  of  obtaining  o'  under  the  surface  for  an  input  angle  of  0.  The  result 
of  this  integration  is: 

'2  ' nl  \2  ”1 


fl  ^ ^ 7T  R 0 


exp 


- Os)  * (*)  ] 


(18) 


where: 


0 = 


r 2 , 2 A2  1l/Z 

k (/?_  1}  + *0  J 


(19) 


3-1 


NADC-78253-30 


I 


I 


I 


-I 

I 


i 

s 


f 

r 

i 


n* 

In  equation  (19),  P = — . The  output  beam  distribution  function  depends  on  the 

initial  beam  divergence  and  the  sea  state.  Equation  (19)  is  shown  in  graphical 

form  in  Figure  3.  It  can  be  seen  in  Figure  3 that  there  is  very  little  difference 

between  the  8 = 0.000  radians  curve  and  the  8 = 0.010  radians  curve  for  wind 
o o 

speeds  greater  than  2 meters  per  second.  Most  laser  divergence  cases  would  fall 
close  to  these  two  curves.  Note  that  in  Figure  3,  for  a low  divergence  laser  beam 
(0. 000  to  0.010  radians),  the  beam  spreading  will  be  28  milliradians  for  a 4- 
meter-per-second  wind  speed  and  40  milliradians  for  an  8-meter-per-second 
wind  speed.  These  values  are  comparable  to  the  values  in  Table  2 for  the  small 
beam  case  when  the  probability  is  80  percent.  The  80-percent  probability  case 
corresponds  very  closely  to  the  output  beam  distribution  function  being  multi- 
plied by  it  is  for  this  reason  that  we  use  the  80-percent  probability  case. 

This  beam  spreading  result  will  be  incorporated  into  an  equation  for  the  irradi- 
ance  in  the  next  section. 
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SECTION  4 
BULK  SCATTERING 


In  this  section  the  approach  of  Arnush  (Reference  5)  is  used  to  obtain  an 
expression  for  the  irradiance  as  a function  of  depth  z.  The  small  angle  scatter- 
i.ig  approximation  of  electron  scattering  theory  is  used.  The  reasonably  good 
fit  at  small  angles  to  existing  experimental  data  is  shown  in  Figure  4.  The 
normalized  scattering  function  2 ( 6 ),  which  is  defined  as  the  volume  scattering 
function  divided  by  the  scattering  coefficient,  is  of  the  following  form: 


SW  = 


y exp  (-  y 9) 
2 *9 


y » 10 . 


(20) 


The  integral  of  the  normalized  scattering  function  is  used  in  a formal  expression 
for  the  ir radiance.  The  result  of  this  calculation  is: 


h (z,  r,  t) 


JL& 


where: 


2 s z' 


3 y 


2 

+ z 


2 

0 . 


(21) 


(22) 


In  the  above  equations, 

h = the  scalar  irradiance, 

g(t)  = a time -dependent  amplitude  which  depends  on  the  source  strength, 
a = the  volume  absorption  coefficient, 
s = ths  volume  scattering  coefficient, 

R0=  the  radius  at  the  surface  at  which  the  beam  decays  to  - times  the 
value  at  the  center,  ?nd 

0 = the  angular  spread  of  the  beam  just  below  the  surface. 
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Figure  4.  Small  Angle  Scattering  Approximation 
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The  amount  of  energy  passing  through  an  aperture  of  radius  p centered  on 
the  beam  is  given  by: 


E (p,  z,  t) 


h (z,  r,  t) 


(23) 


—a  rr 

= g(t)  e 


(24) 


In  equation  (24) , one-half  the  energy  passes  through  the  aperture  when: 
2 


exp 


= 1/2 


R 


(25) 


1 


or: 


1/2 


2 1/2  2 1/2 
= (0.693  Rx)  a 0.833  (Rj  ) 


(26) 


And,  finally, 


,1/2 


Pl/2  ■ °-M3l77-  * Ro  + 22  ®2’ 


(27) 


Equation  (27)  yields  the  functional  relationship  between  the  half -power  radius 
( f*i /2)  311(1  the  beam  dissipation  mechanisms.  The  first  term  in  the  parentheses 

O 

shows  the  z dependence  due  to  bulk  scattering.  This  term  will  clearly  domin- 
ate at  large  depths.  The  second  term  is  the  surface  spot  size  and  dominates  at 
small  z values.  The  last  term  represents  the  surface  spreading  due  to  beam 
divergence  and  sea  state. 

Bv  equating  the  spot  size  term  and  the  spreading  term  in  equation  (27)  to 
the  bulk  scattering  term,  the  depth  at  which  bulk  scattering  begins  to  dominate 
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I 


1 


can  be  calculated.  The  depth  to  which  the  beam  radius  term  dominates 
over  the  scattering  term  is: 


(28) 


For  y = 10,  Rq  = 10  meters,  and  s = 0.1  meters  \ we  find  that  50  meters. 
Therefore,  the  beam  radius  term  dominates  the  scattering  term  to  a depth  of  50 
meters  for  an  initial  beam  radius  of  10  meters.  Clearly,  the  initial  beam  size 
is  very  important  in  determining  the  half -power  radius. 


A comparison  of  the  bulk  scattering  term  and  the  surface  spreading  term 
yields : 


o 


2 


(29) 


where  z^  is  the  depth  at  which  the  bulk  scattering  term  and  the  surface  spread- 
ing term  are  equal.  For  a surface  spreading  term  of  0.040  radians,  s =0.05 
meter  , and  7 ~ 10,  we  find  that  z 0 =*  4. 8 meters.  This  is  a very  important 
result.  The  bulk  scattering  term  dominates  over  the  surface  spreading  term 
from  4. 8 meters  and  deeper,  Qven  for  very  good  water  (s  = 0.05  meter  ) and  an 
8-meter-per-second  wind  speed. 

Various  plots  of  equation  (27)  are  shown  in  Figure  5 for  a scattering  co- 
efficient of  0.1  meter-i.  The  parameters  used  for  RQ  and  0 are  not  physically 
realistic  but  are  presented  for  didactic  reasons.  The  benefit  of  using  a small 
beam  for  increasing  the  irradiance  is  readily  apparent.  It  is  also  clear  that 
there  is  very  little  benefit  in  correcting  for  the  sea  surface;  at  50  meters  the 
correction  would  result  in  only  a 2-percent  decrement  in  beam  radius  for  the 

conditions  given  in  Figure  5.  The  dominance  of  the  bulk  scattering  term 
3/2 

(pa  z ) at  large  depths  is  quite  clear  in  Figure  5. 
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Figure  5.  Beam  Radius  Versus  Depth 
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SECTION  5 

EXPERIMENTAL  APPARATUS 


The  experimental  studies  were  begun  with  the  idea  in  mind  that  adaptive 
optic  techniques  might  be  utilized  to  compensate  for  the  fine-scale  roughness 
of  the  sea  surface.  It  is  the  fine-scale  roughness  due  to  capillary  waves  that 
contributes  markedly  to  beam  spreading.  Therefore,  it  was  decided  (1)  to 
measure  the  capillary  wave  heights,  as  this  defines  the  needed  alteration  of 
the  wave  front;  (2)  to  measure  the  frequency  content  from  which  the  length  of 
the  fine-scale  roughness  and  the  time  constant  of  the  surface  state  can  be  cal- 
culated; and  (3)  to  measure  the  wave  slope  distribution,  whi^H  determines  the 
beam  divergence  under  the  surface. 

The  experimental  apparatus  for  measuring  wind-driven  wave  heights  is 
shown  in  Figure  6.  The  wavetank  is  a rectangular  box  (28  cm  width  x 19  cm 
height  x 117  cm  length)  fabricated  from  clear  Plexiglas.  Mounted  in  one  end 
of  the  wavetank  is  a tangential  blower  driven  by  a variac.  The  blower  forces 
air  into  an  aluminum  plenum  which  provides  a laminar  flow  of  air  over  the 
water  surface  for  a length  of  20  cm.  The  capillary  waves  are  generated  within 
this  region  of  uniform  air  flow.  The  air  velocity  across  the  exit  opening  of  the 
plenum  varies  less  than  10  percent  for  the  range  of  wind  speeds  used  in  these 
experiments  (2.2  meters  per  second  to  7. 6 meters  per  second).  The  wind 
speeds  were  measured  with  a Pitot  tube  and  an  Alror  Velometer.  At  the  end 
of  the  wavetank  opposite  to  the  blower  is  a sloping  aluminum  plate  used  to  damp 
out  the  capillary  waves.  A Spectra  Physics  Model  155  He-Ne  laser  is  mounted 
vertically  above  the  wavetank  and  illuminates  the  water  surface  10  centimeters 
from  the  end  of  the  plenum.  The  laser  beam  diameter  is  1 millimeter.  The 
laser  beam  traverses  the  wavetank  and  enters  the  light  collection  system  through 
a narrow  slit.  The  slit  size  for  the  wave  height  measurements  was  0.6  cm  x 
1.2  cm.  A slit  of  this  size  eliminates  the  light  which  has  traversed  long  paths 
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Figuie  6.  Wav  stank  Apparatus 
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through  the  water  from  being  collected.  The  light  collection  system  includes 
a diffusing  screen  and  a collecting  lens  which  uniformly  illuminates  the  surface 
of  an  RCA  6199  photomultiplier  tube.  This  system  eliminates  the  problem  asso- 
ciated with  the  photomultiplier  tube  output  varying  with  the  location  of  the  illum- 
ination. The  photomultiplier  tube  output  is  amplified  and  simultaneously  dis- 
played by  a Tektronix  Type  547  Oscilloscope  and  is  recorded  by  a Honeywell 
Model  1508  Visicorder.  The  visicorder  galvanometer  has  a 1 kHz  response. 

The  wavetank  was  filled  to  a depth  of  8. 5 centimeters  with  distilled  water.  The 
water  was  filtered  to  eliminate  scattering  from  particulates.  For  the  wave 
height  measurements,  ink  was  added  to  the  water  to  yield  a solution  attenuation 
coefficient  of  « = 0.85  cm-1.  This  coefficient  was  obtained  by  filling  a 1- 
centimeter  sample  cell  with  the  ink  solution  and  measuring  the  attenuation  of 
an  He-Ne  laser  beam  through  the  cell  with  an  EG&G  Radiometer. 

The  apparatus  used  to  measure  the  wave  spectrum  is,  to  a large  extent, 
the  same  as  that  in  Figure  6.  The  differences  are  that  clear  water  is  used,  the 
collection  slit  is  narrowed  to  0.2  cm  x 1.2  cm,  the  amplifier  is  eliminated,  and 
the  visicorder  is  replaced  by  a Tektronix  3L5  Spectrum  Analyzer  in  a Type  564 
Storage  Oscilloscope. 

For  the  wave  slope  measurements  the  slit  dimensions  were  0.2  cm  x 1.2 
cm.  The  photomultiplier  tube  was  translated  in  0, 2-centimeter  increments 
parallel  to  the  wind  direction.  The  signal  was  averaged  using  an  operational 
amplifier  integrating  circuit  with  a 10-second  time  constant.  At  least  2 minutes 
were  allowed  to  elapse  in  each  position  to  accumulate  data  before  the  voltage 
was  read  by  a Hewlett  Packard  Model  400H  vacuum  tube  voltmeter. 


5-3 


NADC-78253-30 


I 

1 

8 SECTION  6 

| EXPERIMENTAL  RESULTS 

f Fluctuations  in  the  laser  intensity  were  recorded  in  real  time  by  the  visi- 

J corder.  The  wave  height  distribution  was  calculated  at  various  wind  speeds 

j from  variations  in  the  laser  beam  intensity  after  traversing  the  wavetank. 

I 

The  governing  equation  is: 
i dl 

Y = - « 8 (30) 

i where: 

dl  = the  incremental  change  in  intensity  I , 

; -1 

« = the  attenuation  coefficient  0. 85  centimeter  , and 

S = the  trough  to  crest  wave  height. 

Using  the  above  attenuation  coefficient,  a 1 -millimeter  high  wave  yields  an  8.5- 
percent  change  in  signal,  which  is  well  above  the  noise.  In  fact,  for  the  wave 
height  experiments  the  high-frequency  noise  (~200  Hz),  which  averages  out,  is 
'~5  percent,  and  the  low  frequency  noise  (~10  Hz)  is  ~1  percent.  This  permits 
wave  heights  ~ 0. 1 millimeter  to  be  measured.  In  Table  3 the  wave  heights  and 
the  standard  deviations  are  listed  for  various  wind  velocities.  The  standard 


TABLE  3.  WIND-DRIVEN  WAVE  HEIGHTS 


WIND  SPEED 
(meters/second) 

WAVE  HEIGHT 
(millimeters) 

2.2 

0.06+0.15 

3.1 

0.32+0.16 

4.0 

1.2  +0.7 

4.9 

1.5  +0.6 

5.8 

2,8  +0.9 

6.7 

3.9  +1.3 

7.6 

8,8  +3.6 

G-l 
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deviations  on  the  two  lowest  data  points  (2.2  m/sec  and  3.1  m/see)  are  a result 
of  the  system  noise.  The  standard  deviations  on  the  rest  of  the  data  are  a re- 
sult of  the  wave  heights  being  broadly  distributed  at  each  wind  speed.  The  wave 
height  distribution  for  a 5.  8-meter-per-second  wind  speed  is  shown  in  Figure  7. 
The  data  in  Figure  7 typify  the  results  at  other  wind  speeds.  Therefore,  each 
wind  speed  is  not  well-characterized  by  a single  wave  height.  A plot  of  the  data 
in  Table  3 is  shown  in  Figure  8.  It  can  be  seen  in  Figure  8 that  the  mean  value 
of  the  wave  height  increases  monotonically  with  wind  speed  from  a value  of 
S =0.06  mm  at  2.2  meters  per  second  to  I = 8.8  mm  at  7.6  meters  per  second. 
Several  features  in  Figure  8 are  worth  noting.  The  highest  mean  wave  height 
(S  = 8. 8 mm)  is  within  the  allowable  height  for  capillary  waves  (1.73  cm  x 
0.73  = 1.24  cm).  The  second  feature  is  that  capillary  waves  first  appear  at  a 
wind  speed  of  2.2  meters  per  second.  This  value  is  slightly  on  the  low  side  for 
our  fetch  (30  cm)  if  the  data  are  compared  with  Cox's  (Reference  6).  Cox  found 
that,  for  a 1.47 -meter  fetch,  UQ  =2.9  meters  per  second,  where  UQ  is  the 
wind  speed  below  which  the  mean  square  slope  is  extremely  small.  In  the 
present  study  the  extrapolation  of  the  mean  square  slope  to  zero  yields  a higher 
value  for  the  onset  wind  speed  for  capillary  wave  generation  as  will  be  seen 
subsequently. 

The  spectral  intensity  versus  frequency  curves  for  specific  wind  speeds 
are  shown  in  Figure  9.  It  can  be  seen  in  these  plots  that  the  lower  frequencies 
fall  off  in  intensity  with  increasing  wind  speed,  while  the  higher  frequencies 
increase  in  intensity.  Most  of  the  spectral  intensity  occurs  below  100  Hz  (A> 
0.35  centimeter)  even  at  the  highest  wind  speeds  measured  in  these  experiments 
(7.6  meters  per  second) . There  are  very  few  waves  of  a frequency  greater  than 
300  Hz;  this  corresponds  to  wavelengths  less  than  0.17  centimeters.  Therefore, 
in  these  experiments  most  of  the  measurements  are  on  waves  in  the  wavelength 
range  of  0.35  centimeters  to  1.7  centimeters  and  in  the  frequency  range  100  Hz 
to  13  Hz. 
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Figure  9.  Wave  Spectral  Content  at  Specific  Wind  Speeds 
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In  Figures  10  through  16,  the  wave  slope  distributions  are  shown  at  wind 
speeds  in  the  range  2.2  meters  per  second  to  7.6  meters  per  second.  The  data 
in  each  figure  are  corrected  for  photomultiplier  tube  dark  noise  and  ambient 
illumination.  In  monitoring  the  pulse  heights  during  the  integration  time,  it 
was  found  that  there  was  a systematic  but  small  fall-off  of  signal  level  with 
angle;  the  data  as  shown  are  also  corrected  for  this.  The  slit  width  (2  milli- 
meters) corresponds  to  1.22  degrees  of  collection  angle  which  is  equal  to  ~4. 9 
degrees  of  wave  slope.  The  tube  is  translated  in  increments  of  2 millimeters; 
therefore,  there  is  no  overlap  in  data  collecting  from  one  reading  to  the  next. 

In  Figure  10  the  histogram  for  the  2.2  meters  per  second  wind  speed  is 
shown.  There  are  waves  present  at  this  wind  speed  as  indicated  from  the  wave 
height  measurements.  However,  all  of  the  slopes  present  are  within  the  central 
4.9  degrees  (+  2.45  degrees);  therefore,  no  intensity  is  measured  beyond  this 
region.  It  would  take  a smaller  slit  to  measure  the  slope  distribution  at  this 
wind  speed  (2. 2 meters  per  second).  An  additional  feature  of  the  data  is  that 
the  areas  under  the  histograms  in  Figures  10  through  16  are  nearly  constant. 
(Note  the  scale  change  in  Figures  12  through  16. ) The  nearly  constant  areas 
result  from  the  fact  that  nearly  all  of  the  refracted  light  is  collected.  Therefore, 
at  higher  wind  speeds  the  central  region  diminishes  as  the  intensity  at  higher 
angles  grows.  In  looking  at  Figures  10  to  16  it  is  noticed  that,  for  wind  speeds 
equal  to  4. 9 meters  per  second  and  greater,  the  peak  of  the  slope  distribution 
is  shifted  towards  the  downwind  direction.  This  same  behavior  was  seen  by 
Cox  and  Munk  (Reference  4)  and  Prettyman  and  Cermak  (Reference  7)  in  their 
ocean  measurements.  Besides  this  asymmetry  in  the  data,  it  is  seen  that  the 
intensities  are  different  at  high  angles  (greater  than  30  degrees)  for  the  wind- 
ward and  leeward  sides  of  the  waves.  The  shapes  of  the  histograms  show  that 
the  shape  of  the  wind-driven  waves  is  affected  by  the  wind.  On  the  windward 
side  of  the  waves,  the  wave  slope  is  in  the  low  angle  positions  longer  and  in 
the  high  angle  positions  for  a shorter  time  than  on  the  leeward  side  of  the  waves. 
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From  the  data  in  Figures  10  to  16,  the  rms  wave  slope  may  be  calculated 
for  each  wind  speed.  A plot  of  this  calculation  is  shown  in  Figure  17.  The 
rms  wave  slope  rises  sharply  with  wind  speed  up  to  5 meters  per  second,  then 
begins  to  level  off.  If  the  steeply  rising  portion  of  the  wave  slope  versus  wind 
speed  curve  is  extrapolated  to  zero,  an  onset  velocity  is  found  of  UQ  ~ 3 meters 
per  second.  This  value  (U0  ~3  meters  per  second)  for  the  short  fetch  used  in 
this  study  is  quite  compatible  with  the  previous  numbers  given  for  the  Cox  study 
(Reference  6).  In  Figure  17  the  data  of  Cox  are  plotted  for  a fetch  of  1.47  meters 
for  comparison.  The  values  found  in  the  present  study  for  the  rms  wave  slopes 
rise  more  steeply  with  wind  speed  and  level  off  at  a higher  value.  The  ocean  data 
of  Cox  and  Munk  (Reference  4)  obtained  from  sunglint  measurements  are  also  shown 
in  Figure  17.  The  ocean  data  differs  from  the  wavetank  data  in  two  major  re- 
spects: (1)  there  is  no  onset  velocity  in  the  2 to  3 meters  per  second  range,  and 
(2)  the  rms  values  for  the  wave  slopes  at  a given  wind-speed  are  lower  in  the 
open  ocean  than  in  wavetanks.  The  smaller  wave  slopes  in  the  open  ocean  have 
been  attributed,  at  least  in  part,  to  the  suppression  of  capillary  wave  growth  by 
gravity  waves  (Reference  8). 

In  Figure  18  the  full -angle  beam  divergence,  which  contains  50  percent  of 
the  energy  of  the  beam,  is  plotted  against  wind  speed  for  three  different  cases. 

The  first  case  is  represented  by  the  four  data  points  of  Prettyman  and  Cermak 
(Reference  7)  which  show  the  beam  divergence  to  be  about  4 degrees  and  inde- 
pendent of  wind  speed.  Their  data  was  taken  in  the  ocean  using  a laser  source 
illuminating  a submerged  screen.  The  location  of  the  laser  spot  (l/8-inch 
diameter)  was  filmed.  The  data  from  the  present  study  are  also  plotted  in 
Figure  18;  the  data  points  were  obtained  from  Figures  10  through  16  by  finding 
the  full  width  of  each  histogram  which  contains  half  the  power.  The  theoretical 
curve  is  calculated  from  a one-dimensional  model  using  the  value  for  the  rms 
slope  from  the  Cox  and  Munk  study  (Reference  4).  There  is  a fair  amount  of 
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Figure  17.  RMS  Wave  Slope  Versus  Wind  Speed 
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Figure  18.  Beam  Divergence  Versus  Wind  Speed 
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discrepancy  in  the  three  results.  The  wavetank  measurements  cf  the  present 
study  are  not  expected  to  represent  ocean  conditions  because  of  the  short  fetch 
and  the  absence  of  gravity  waves.  However,  it  is  expected  that  there  would  be 
better  agreement  between  the  laser  data  of  Prettyman  and  Cermak  and  the  theo- 
retical model  which  incorporates  the  Cox  and  Munk  data.  Also,  a stronger  de- 
pendence of  laser  divergence  on  wind  speed  is  expected  from  the  present  study 
than  is  shown  by  the  Prettyman  and  Cermak  data. 

One  of  the  reasons  for  the  lower  rms  wave  slopes  in  the  Cox  and  Munk  paper 
could  be  related  to  the  geometry  of  the  capillary  waves.  In  Figure  19  the  capil- 
lary waveforms  as  calculated  by  Crap  per  (Reference  3)  are  shown  for  a/x 
(amplitude/wavelength)  values  equal  to  0.73,  0.53,  and  0,34.  When  considering 
reflecting  light  from  the  a/x  = 0.73  capillary  waves,  it  is  found  that  adjacent 
waves  obscure  the  higher  slopes  from  being  measured.  For  example,  for  the 
case  a/x  = 0.73  and  the  elevation  angle  of  the  light  source  is  90  degrees  (directly 
overhead),  the  highest  measurable  slope  angle  is  approximately  40  degrees.  If 
the  above  wave  is  assumed  to  be  planar,  then  wave  slopes  greater  than  or  equal 
to  40  degrees  account  for  34  percent  of  the  reflection  and  refraction.  This  num- 
ber is  calculated  by  taking  the  horizontal  projection  of  the  wave  envelope  and 
measuring  the  fraction  (F)  of  the  projection  for  which  the  slope  angle  (tr)  is 
greater  than  or  equal  to  40  degrees.  For  a/x  = 0.73,  F = 51  percent  for  <?  } 

30  degrees.  For  a/x  = 0.53,  F = 39  percent  for  <j  ^,30  degrees.  And,  for  a/x  = 
0.34,  F = 32  percent  for  * ^,30  degrees.  In  this  paragraph,  extreme  examples 
(high  a \ ratios)  have  been  taken  to  make  the  point  that  the  reflection  technique 
systematically  prevents  the  higher  slope  angles  from  being  measured. 

As  seen  in  Figure  8,  the  mean  capillary  wave  height  for  the  highest  wind 
speed  (7.6  meters  per  second)  is  8.8  millimeters;  therefore,  for  capillary  waves 
( x ^ 1.7  centimeters)  the  ratio  a/X  ^0.52,  Consequently,  a high  percentage 
(F  ~ 39  percent)  of  wave  slopes  greater  than  30  degrees  should  be  observed  in 
Figure  16. 
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However,  the  percentage  measured  from  Figure  16  is  F«8  percent.  This  dis- 
crepancy could  be  due  to  impurities  in  the  water  used,  the  effects  of  wind  on  the 
wave  shape,  or  applying  a two-dimensional  model  to  a three-dimensional  world. 

In  comparing  these  data  to  the  model  by  Crapper  (Reference  3),  there  is  one 
last  point  worth  making.  The  highest  wave  height  measured  in  all  of  our  data 
taking  was  1.23  centimeters;  if  it  is  assumed  that  this  wave  height  was  associated 
with  a capillary  wave  of  maximum  length  (1.7  centimeters),  then  a/X  = 0.72. 

This  value  is  in  excellent  agreement  with  the  theoretical  value  (a/X  = 0. 73)  of 
Crapper. 
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SECTION  7 
CONCLUSIONS 


The  following  conclusive  can  be  drawn  from  our  study: 

a.  The  beam  spreading  due  fo  the  roughness  of  the  air-sea  interface  is 
small  compared  with  the  spreading  due  to  bulk  scattering.  Therefore, 
compensating  for  the  interface  would  be  of  little  benefit. 

b.  Using  a small-diameter  laser  beam  substantially  increases  the  irra- 
diance  on  target  for  depths  of  50  meters.  and  less. 

c.  The  mean  heights  for  capillary  waves  are  in  the  range  0.06  to  8.8 
millimeters  for  wind  speeds  2.2  to  7.6  meters  per  second.  These 
heights  are  consistent  with  theoretically  allowed  values  for  capillary 
waves. 

d.  The  spectral  content  of  the  wind-driven  capillary  waves  is  in  the  range 

0 to  400  Hz  with  most  of  the  intensity  below  100  Hz.  The  high -frequency 
portions  of  the  spectral  curves  grow  at  the  expense  of  the  low-frequency 
portions  as  the  wind  speed  increases. 

e.  The  wavetank  data  on  wave  slopes  yield  higher  values  for  full-width 
beam  divergence  («-5  degrees)  than  open-ocean  data  (~  2 degrees  and 
~ 4 degrees)  at  wind  speeds  in  the  range  of  5 to  8 meters  per  second. 

f.  Fewer  wave  slopes  of  high  angle  are  measured  in  these  experiments 
than  are  predicted  by  theory, 

g.  The  onset  wind  velocity  for  wave  generation  is  found  to  be  ~ 3 meters 
per  second. 

h.  The  shape  of  the  wind-driven  waves  is  affected  by  the  wind, 

i.  The  wave  height  measurements  furnish  experimental  support  to  the 
theoretical  model  of  Crapper  (Reference  3). 
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